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ABSTRACT 

PSR J1738+0333 is one of the four millisecond pulsars known to be orbited by a white 
dwarf companion bright enough for optical spectroscopy. Of these, it has the shortest 
orbital period, making it especially interesting for a range of astrophysical and gravity 
related questions. We present a spectroscopic and photometric study of the white dwarf 
companion and infer its radial velocity curve, effective temperature, surface gravity 
and luminosity. We find that the white dwarf has properties consistent with those 
of low-mass white dwarfs with thick hydrogen envelopes, and use the corresponding 
mass-radius relation to infer its mass; M WD = 0.18]+o'.oo5 M©. Combined with the 
mass ratio q = 8.1 ± 0.2 inferred from the radial velocities and the precise pulsar 
timing ephemeris, the neutron star mass is constrained to Mpsr = 1.47^g'og Mq. 
Contrary to expectations, the latter is only slightly above the Chandrasekhar limit. 
We find that, even if the birth mass of the neutron star was only 1.20 M Q , more than 
60% of the matter that left the surface of the white dwarf progenitor escaped the 
system. The accurate determination of the component masses transforms this system 
in a laboratory for fundamental physics by constraining the orbital decay predicted by 
general relativity. Currently, the agreement is within lex of the observed decay. Further 
radio timing observations will allow precise tests of white dwarf models, assuming the 
validity of general relativity. 

Key words: binaries: close - pulsars: general - stars: neutron - white dwarfs - 
processes: gravitation - individual: PSR J1738+0333 



1 INTRODUCTION 

Millisecond pulsars (MSPs) are extreme in many ways. Their 
interior consists of the densest form of observable matter 
known and they can spin at least as fast as 716 times per 
second ( Hessels et al.|[2~006 |. Hence, they offer a rich labo- 
ratory for a wide range of physical inquiry: Mass measure- 
ments provide direct comparison to quantum chromodynam- 



field tests of general relativity (Taylor & Weisberg 1982 



ics' predictions for the state of ultra-dense matter ( Lattimer 
|fc Prakash||2004| |Demorest et al.||2010| ) and studies of their 
orbits in binaries have provided the first confirmation for 
gravitational wave emission and the most stringent strong- 
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Weisberg et al.|2010| [Kramer et al.|2006 l 

Most of the fastest spinning Galactic-disk pulsars are 
paired with low mass helium-core WDs (hereafter LMWDs, 
for recent reviews see |Lorimer|2008} |Tauris|2?JlT| , and their 
fast spins and weak magnetic fields are thought to be the 
product of mass transfer from the progenitor of the WD, 
a process also known as recycling. As the progenitor star 
evolves, it fills its Roche lobe and loses its envelope, either 
while on the main sequence (for sufficiently short initial pe- 
riods), or when moving up the red-giant track (Webbink 



et al. 19831. The mass transfer rate is a strongly increas- 



ing function of the initial orbital period and donor mass 
( Tauris fc Savonije|1999 1, and is expected to be at a stable, 
sub-Eddington rate (< 10 -8 Mq yr _1 ) for light companions 
in relatively tight orbits. The final result of such long-term 
(nuclear timescale) mass transfer is a highly circular (due 
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to fast tidal dissipation in the secondary) close binary con- 
sisting of a fast spinning MSP and a low mass, helium-core 
WD. 

These systems are important for several reasons. First, 
it is these binaries that allow one to probe certain aspects 
of the radiative properties of gravity that are poorly con- 
strained by the relativistic effects seen in double neutron 
stars, like the Hulse-Taylor or the double pulsar. For exam- 
ple, in a wide range of theories, the rate of gravitational wave 
emission is driven by a leading dipolar term that depends 
crucially on the difference in gravitational binding energies 
between the binary members. Hence, if accurate component 
masses can be determined, one can directly confront the 
predictions of different gravity theories in terms of dipolar 
radiation with observations. 

Second, measuring their masses provides access to the 
accretion process and evolution of these systems as well as 
the formation of MSPs, the only neutron stars with secure 
precise masses significantly above the Chandrasekhar limit 
( Freire et al.|2011 Demorest et al.| 20f0). In addition, obser- 
vational constraints on the upper mass limit of stable neu- 
tron stars, constrains the equation of state for super-dense 
matter. 

Unfortunately, precise MSP and companion masses can 
be determined from timing in exceptional cases only: either 
when the orbit is (unexpectedly) eccentric, allowing for a 



measurement of the rate of advance of periastron (Freire 
et al.|20rTj), or if the system has an orbit seen edge on fKaspi 



et al.|1994| |Jacoby et aL|2005| |Demorest et al.|2010[ ) which 
allows for a measurement of pulse time-of-arrival (TOA) de 
lays due to the curvature of space-time around the compan- 
ion (Shapiro delay, Shapiro 1964 1 . 

Fortunately, another method exists that relies on com- 



bined optical and radio timing observations (van Kerkwijk 
|et al.||E)96l |Callanan et~aL]pt8l ). If the WD companion is 
bright enough for detailed spectroscopy, a comparison of its 
spectrum with model atmospheres yields its effective tem- 
perature and surface gravity. These can then be compared 
with a mass-radius relation for LMWDs to yield its mass. 
Combining the radial velocity for the white dwarf with the 
pulsar timing measurements yields the mass ratio and there- 
fore the mass of the pulsar. 

In a companion paper (Van Kerkwijk et al. 2011, in 
prep.; vK+12 hereafter), we test this method on PSR 
J1909— 3744, for which the masses are precisely known from 
timing. We find it reliable and are confident to apply it also 
to other similar systems. Here, we report on an application 
to PSRJ1738+0333, a pulsar-LMWD binary in a 8.5 hours 
orbit. This system is particularly interesting because of its 
short period, which places it in a poorly studied regime for 
accretion physics of neutron stars, where nuclear-driven evo- 
lution competes with magnetic braking and gravitational 



radiation (c.f. Phinney & Kulkarni 1994). Furthermore, the 



short orbital period implies relatively rapid orbital decay, 
making the system interesting for radiative tests of gravity, 
for which, as described above, knowledge of the component 
masses is necessary. 

The paper is organized as follows: In §[2] we report re- 
sults from radio timing, necessary for calculations through- 
out the rest of the paper. These are presented in detail sep- 
arately ( jFreire et al.|[20T2} from now on Paper II). In §[3] 
we describe the spectroscopic and photometric observations 



and in §|4] we present our results. We discuss our findings 
and comment on the evolution of the system and its impor- 
tance for gravity tests in §[5] Finally, in §[6] we summarize 
our results. 



2 PSRJ1738+0333: RADIO OBSERVATIONS 

PSR J 1738+0333 was discovered in a 20-cm high Galactic 
latitude survey in 2001 ( Jacoby 2005]), carried out with the 
multi-beam receiver of the Parkes Telescope. The pulsar has 
a spin period of 5.85 ms and orbits a low-mass helium-core 
WD companion in a 8.5 h orbit. Since 2003 it has been reg- 
ularly timed with the 305 m Arecibo Telescope, leading to 
~ 17000 times of arrival with typically 3 /its uncertainties. 
The corresponding timing solution provides measurements 
of the system's parallax and proper motion, and a signifi- 
cant detection of the intrinsic orbital period derivative (see 
Paper II for details). In Table[2] we list the measured spin, 
Keplerian and astrometric parameters of the system. 

The spin period derivative is that of a typical low- 
surface magnetic field pulsar (Bo = 3.7 x 10 8 G), and the 
characteristic age (= P/2P) after subtracting the kine- 
matic effects (Paper II) is 4.1 Gyr. The parallax measure- 
ment corresponds to a distance of d — 1.47 ± 0.10 kpc. 
The system's proper motion combined with the parallax 
implies transverse velocities of v a — dfi a = 49kms _1 and 
vg — dfis = 36kms _1 in a and S respectively. In sjZ] we 
combine these values with the systemic radial velocity, 7, 
to derive the 3D spatial velocity and calculate the Galactic 
orbit of the binary. The estimate for the orbital eccentricity 
is one of the lowest observed in any binary system: When 
Shapiro delay is accounted for in the solution (Paper II), the 



apparent eccentricity diminishes to e = (3.5 ± 1.1) x 10 . 
We discuss the implication of this for evolutionary scenarios 
in £0 



3 OPTICAL OBSERVATIONS 
3.1 Spectroscopy 

Our main data set consists of eighteen long-slit phase re- 
solved spectra of PSR J1738+0333, obtained with the Gem- 
ini South telescope at Cerro Pachon on ten different nights 
between April and June 2006 (see Tablejl]). For our ob- 
servations we used the Gemini Multi-Object Spectrograph 
(GMOS-S). The CMOS detector consists of three 2048 x 
4608 EEV CCDs, each of which was read-out at 2 x 2 binning 
by a different amplifier, giving a scale of 0'.'14 per binned 
pixel in the spatial direction, and, with the 1200 lines per mm 
B1200 grism, 0.4 A per binned pixel in the dispersion direc- 
tion. We chose a relatively wide, l'.'5 slit, to minimize atmo- 
spheric dispersion losses (see below). This meant that the 
resolution was set by the seeing, at ~ 3 A, or ~200kms _1 
at 4300 A. In order to cover the higher Balmer lines, we cen- 
tred the grating at 4300 A, for a wavelength coverage from 
3500 to 5100 A. 

All exposures had integration times of 3720 s and were 
followed by an internal flat-field exposure and a Copper- 
Argon (CuAr) exposure for wavelength calibration. The slit 
was oriented to include a bright comparison star located 
25'.'2 at position angle 127?57 (north through east) of the 
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Table 1. Log of observations and radial velocity measurements. 



fR %D Aw 

Date MJD midbar <j> (kms -1 ) (kms -1 ) (kms -1 ) AB' 

(1) ' (2) (3) ' (4) ' (5) (6) 



Gemini,GMOS-S 



2006 Apr 27 


53852 


310219 





0250 


+56 3 + 


g 


— 209 ± 27 


— 265 + 27 


2 


91 


± 





06 




53852 


366314 





1831 


+49.1 ± 


g 


— 143 ± 26 


— 192 +_ 26 


2 


88 


± 


o 


06 


2006 Apr 28 


53853 


295453 





8019 


+53.7 ±0 


5 


-100 ± 14 


-154 ± 14 


2 


88 


± 





04 




53853 


350638 





9575 


+68.8 ±0 


6 


-185 ± 15 


-254 ± 15 


2 


88 


± 





04 


2006 May 07 


53862 


333037 





2749 


+40.1 ± 


5 


-35 ± 13 


-75 ± 13 


2 


83 


± 





03 




53862 


391933 





4409 


+84.9 ±0 


6 


+162 ± 21 


+77 ± 21 


2 


87 


± 





01 


2006 May 26 


53881 


198674 





4489 


+60.8 ± 1 


1 


+121 ± 37 


+60 ± 37 


3 


03 


± 





11 




53881 


252933 





6018 


+67.0 ±0 


9 


+33 ± 32 


-34 ± 32 


2 


97 


± 





09 


2006 May 27 


53882 


352291 





7005 


+43.1 ± 


5 


-6 ±15 


-49 ± 15 


2 


85 


± 





04 


2006 May 28 


53883 


296760 





3625 


+54.9 ±0 


6 


+28 ± 16 


-27 ± 16 


2 


86 


± 





03 




53883 


350144 





5130 


+53.9 ±0 


6 


+134 ± 17 


-189 ± 17 


2 


85 


± 





03 


2006 Jun 19 


53905 


174549 





0264 


+41.5 ± 


5 


-226 ± 12 


+80 ± 12 


2 


85 


± 





04 


2006 Jun 23 


53909 


170100 





2881 


+39.9 ±0 


7 


-5 ±14 


-45 ± 14 


2 


88 


± 





03 




53909 


210618 





4023 


+95.6 ± 1 


1 


-103 ± 36 


-45 ± 36 


2 


84 


± 





04 


2006 Jun 26 


53912 


156147 





7045 


+60.4 ±0 


5 


-7± 14 


-67 ± 14 


2 


85 


± 





04 




53912 


209838 





8558 


+60.7 ±0 


6 


-136 ± 15 


-197 ± 15 


2 


88 


± 





04 


2006 Jun 27 


53913 


120000 





4212 


+42.5 ± 


5 


+79 ± 12 


+37 ± 12 


2 


84 


± 





04 




53913 


176660 





5809 


+43.6 ± 


5 


+106 ± 13 


+62 ± 13 


2 


82 


± 





04 


Keck, LRIS 




























2008 Aug 04 


54682 


377697 





6224 


50 ± 1,+61 ± 5 


-2 ±9 


-52 ±9 


3 


01 


± 





05 



Notes: (1) refers to the barycentric mid-exposure time. (2) is the orbital phase using the ephemeris in Table[2] (3) is the comparison's 
velocity in respect to the solar system barycenter and (4) the raw barycentric velocities of PSR J1738+0333. (5) is the differential 
velocity used to determine the orbit in ^4] Finally, (6) are the differential spectrophotometric magnitudes in B' (equal to B, but limited 
to the wavelength range covered by our spectra; see §31. Here, the errors are the quadratic sum of the photometric uncertainties of the 
WD and the comparison. For LRIS, two velocities are listed for the comparison star, for the blue and red arm, respectively. For the 

white dwarf, the velocity is for the blue arm (see text). 



WD (see Fig. [TJ. We use this star as a local velocity and 
flux standard (since GMOS-S does not have an atmospheric 
dispersion corrector, slit losses vary with offset from the par- 
allactic angle). 

The conditions during the observations were mostly 
good to photometric, but some exposures were taken 
through thin cirrus. The seeing ranged from 0'.'6 to l'/2. For 
flux calibration, we acquired additional frames of the com- 
parison star and the spectro-photometric standard EG 274 
through a 5'.'0 slit on the night of 2006 April 27 (which was 
photometric and had 0'.'8 seeing). Furthermore, for absolute 
velocity calibration, we observed the radial velocity standard 
WD 1743-132 on 2006 June 19. 

The data were reduced using standard and custom rou- 
tines inside the Munich Image and Data Analysis System 
(midas). First, the bias level of each exposure was removed 
using average values from the overscan region. Subsequently, 
we corrected the raw counts on the red and middle chips for 
the small, few percent variations in gain (see vK+12 for 
details on the method), that affected several sets of expo- 
sures (but fortunately not those of the night the flux cal- 
ibrator was taken). Finally, the frames were corrected for 
small-scale sensitivity variations using normalised lamp ex- 
posures, where the normalisation was done both along each 
wavelength position as well as along each spatial position. 
These normalisation steps were required since the lamp spec- 
tra showed rather sharp bumps in the dispersion direction 
whose position and shape was different from bumps seen in 
target spectra, and also varied between sets of spectra (pos- 



sibly because the holographic grating was not illuminated 
exactly identically between the different exposures), while 
in the spatial direction they showed striations due to irreg- 
ularities in the slit. 

For sky subtraction, we selected a 100" region centred 
on the WD, but excluding 5" spots around it and the com- 
parison star. Each column in the spatial direction was fitted 
with a second degree polynomial and the interpolated sky 
contributions at the positions of the WD and the comparison 
were removed. 

Optimally weighted spectra and their uncertainties were 



extracted using a method similar to that of Home ( 1986 1 



The extraction was done separately in each chip and the 
spectra were merged after flux calibration. 

The dispersion solution was established using the CuAr 
spectra taken after each exposure. First, the ID lamp spec- 
trum was extracted by averaging the signal over the spatial 
direction in areas of the chip that coincided with each star. 
Then the lines' positions were measured and identified and 
the dispersion relation was approximated with a 3 rd degree 
polynomial that gave root-mean-square residuals of less than 
0.04 A for typically 18 lines. 

The wide-slit spectra of EG 274 and the comparison 
star were extracted with the same procedure and used to 
calibrate the narrow-slit exposures. Initially, all wide and 
narrow slit data were corrected for atmospheric extinction 
using the average extinction table for La Silla (which should 
be a good approximation to that of Cerro Pachon). Then, 
we calculated the wavelength-dependent flux losses due to 
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5 




Right ascension (J20Q0) 



2048 x 2048 pixels 24 fim on the side, which we read-out 
unbinned. 

The spectra were extracted and calibrated as above. 
Here, on the blue arm we replaced the poorly exposed part 
of the fiat fields shortward of 4000 A with unity and nor- 
malized the rest using a third degree polynomial. On the 
red side we normalized the flat field using a bi- linear fit. 
Wavelength calibration was done using arc spectra and sky 
lines. On the blue arm we used the well exposed arc frames 
taken at the beginning of the night to establish an overall 
solution that had rms residuals of 0.16 A for 22 lines fitted 
with a third-degree polynomial and then calculated offsets 
using the less well-exposed arc frames taken throughout the 
night. For the red arm we used the well exposed arc-frames 
taken interspersed with the science exposures. Here, we cor- 
rected for offsets by shifting the bright OH and O2 lines 
at 8344.602, 8430.174 and 8827.096 A to laboratory values. 
Flux calibration was again done as above; we found that the 
solution was consistent with that obtained from Gemini (see 
also below). 



Figure 1. Finding chart for PSR J1738+0333 (using the SOAR 
V image). Indicated are the white-dwarf counterpart, the slit ori- 
entation used, and the comparison star that was included in the 
slit. 



the finite size of the slit by comparing the wide-slit spec- 
trum of the comparison with each of the narrow-slit spectra. 
The relation was analytically approximated with a quadratic 
function of wavelength that was then applied to the narrow 
slit observations. Finally, the CMOS instrumental response 
was calculated by dividing the spectrum of EG 274 with 
a synthetic template and smoothly interpolating the ratio. 
The template was created by normalizing an appropriate 
DA model atmosphere to the catalogued flux (V = 11.03, 
|Zwitter et aLl[20(Ml see 84.5 for more details on the model 
atmospheres used in this work). Prior to comparison, we 
smoothed the template with a Gaussian kernel to match the 
resolution of the observed spectrum and excluded the cores 
of the Balmer lines. 

Given the possible issues with the detector gain and 
the flat fielding, both of which could affect the flux cali- 
bration, we obtained an additional smaller set of spectra of 
the WD companion, the comparison star and the spectro- 
photometric standard FeigellO using the two-armed Low 
Resolution Imaging Spectrometer (LRIS Oke et al.|1995[ ) of 
the Keck telescope on the night of 2008 August 3 (Table [l]| . 
During the night the sky was photometric and the seeing 
was ~0'.'8. 

For the observations we used the atmospheric dispersion 
corrector and both narrow, 0'.'7, and wide, 8'.' 7 slits. The light 
was split with a dichroic at 6800 A and directed on the two 
arms of LRIS (blue and red arm hereafter). On the blue 
arm we used a 600 lines mm -1 grism, blazed at 4000 A, that 
covers 3100-5600 A with a resolution of AA = 3.2 A or Av = 
220 km s . On the red arm we used the 1200 lines mm -1 
grating, blazed at 8000 A, that covers 7600-8900 A at AA = 
2.1 A or Av = 75kms _1 . The blue-side detector is a mosaic 
of two Marconi CCDs with 4096 x 4096 pixels 15 /an on the 
side, which we read out binned by two in the dispersion 
direction. The red-side detector is a Tektronic CCD with 



3.2 Photometry 

On the night of 2008 February 28, images of the field 
containing PSR J1738+0333 were acquired for us with the 
4.1m Southern Astrophysical Research Telescope (SOAR) 
at Cerro Pachon, Chile, using the Goodman High Through- 
put Spectrograph (Clemens et al.|[2004[ ), with its Fairchild 
4096 x 4096 CCD and B and V filters (with throughputs on 
the Kron-Cousins photometric system). The instrument has 
a plate scale of 0'.'15pix _1 and a usable field of view of 5f0. 
During the run, the sky was photometric and the seeing as 
determined from the images was ~ l'.'8. Two 300 s images 
each in V and B were obtained. Of these, however, the first 
had reduced count rates for all stars and a distorted point- 
spread function, possibly because the telescope and instru- 
ment had not yet settled when the exposure was started; we 
have not used that image. For calibration, sets of 30 s B and 
V images of the photometric standard field PG 1633+099 
were acquired both before and after the science frames. 

Following standard prescriptions, individual frames 
were bias-corrected and flat-fielded using twilight flats. Hot 
pixels and cosmic rays were replaced by a median over their 
neighbours. The instrumental fluxes were measured inside 
3'.'6 radii and then corrected to a radius of 7" using mea- 
surements of bright isolated stars. For the calibration, we 
used 5 standard stars with a range of B — V colors in the 
PG 1633+099 field (Stetson! 119901). Measured magnitudes 



were compared to their catalogued counterparts to derive 
zero-points and colour terms. Both calibration sets yielded 
similar results. Small differences in airmass were corrected 
using standard values for La Silla. The root-mean-square 
residuals of the zero points in both bands were ~0.01 mag. 
We find that the optical counterpart of the WD has V = 
21.30(5) and B = [21.70(7), 21.73(7)], where the two mea- 
surements in B are for the two exposures, and where for the 
errors, we combined in quadrature the measurement and 
zero-point uncertainties. For the comparison star, we mea- 
sure V = 18.00(1) and B = [18.73(2), 18.75(3)]. Since the B 
magnitudes are consistent, we use the averages below. 

We verified our calibration in several ways. First, we 
integrated our flux-calibrated spectra over the B-band fil- 
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ter curve of Bessell (19901. For the comparison star, using 
the wide-slit spectra, we find B' = 18.81 for the Gemini 
spectrum and B = 18.71 for the Keck spectrum. For the 
white dwarf, we find B' = 21.69 for the averaged Gemini 
white dwarf spectra, and B — 21.70 for the single narrow- 
slit Keck spectra. Here, we label the Gemini magnitudes as 
B' , since the GMOS spectra do not fully cover the Bessell-B 
bandpass, which will introduce color terms. 

Second, we tried to calibrate the (/-band GMOS acqui- 
sition images, by calibrating relative to our velocity stan- 
dard, WD 1743-132, which has V = 14.290, B — V = 0.300 



( Mermilliod et al. 19901, and thus, using the relations of 
Fukugita et al.| ( |1996| ), g" = V + 0.56(B-V) -0.12 = 14.34. 
We find g = 18.23 for the comparison star and, using the 
average magnitude difference Ag' = 3.091(17) between the 
WD and the comparison, we infer g' = 21.32 for the WD 
(here, the uncertainty will be dominated by systematics, but 
should be < 0.05 mag). These numbers are consistent with 
the g' = 18.39(7) and 21.42(7) expected from our SOAR 
photometry. 

Looking at individual acquisition frames, the scatter of 
the magnitude difference was ~ 0.05 mag, somewhat larger 
than expected based on measurement noise, though with 
no obvious correlation with orbital phase. We find some- 
what smaller scatter from convolving individual flux cal- 
ibrated WD and comparison spectra with the Bessell B- 
band, and using those to determine differences (see Ta- 
ble[T]). Ignoring the two points from our worst night (2006 
May 26), the root-mean-square scatter is 0.032 mag. Since 
no obvious phase dependence is found, this places a limit 
on the irradiation of the WD atmosphere from the pulsar. 
However, the limit is too weak to be useful: Assuming a 
spin-down luminosity of Lpsr = dE/dt = —Aty 2 I P / P 3 ~ 
4.8 x 10 33 ergss _1 and defining an irradiation temperature 
T irr = (ip S R/47ra 2 cr) 1/4 ~ 3800 K (where from Table [5J 
we inferred a ~ 1.8 x 10 11 cm), the expected orbital mod- 
ulation is only AL/L ~ [TrR 2 VD (LpsR/iTva 2 )/Lwu]smi ~ 
Pta/4Tw D ] sin i < 4 x 10" 3 . 



4 RESULTS 

4.1 Radial velocities 

Radial velocities of the WD, the comparison and the velocity 
standard were extracted by fitting their spectra with tem- 



plates using the method discussed in Bassa et al. ( 2006 1. For 



the comparison, we first classified it using the on-line atlas 
by R. O. GrajQ We find that its spectrum resembles that of 
a GOV star, with an uncertainty of about 1 subtype. Com- 
paring with various spec tra from the UVESPO F|| library of 
high resolution spectra (Bagnulo et al. 20031, we find the 



best fit for the G1V star HD 20807 (where, to match the 
resolution of the observations, we convolve the UVESPOP 
spectra with a Gaussian with FWHM equal to that of the 
seeing, truncated at the slit width). We fitted this template 
to each spectrum for a range of velocities, from —600 to 
600kms _1 with a step size of 5kms _1 . We corrected for 



1 http:/ /nedwww. ipac.caltech.edu/level5/Gray/Graycontent 

2 http:/ /www. sc. eso.org/santiago/uvespop/DATA 



the 11.5kms _1 barycentric velocity of HD 20807 after the 
fact. 

Similarly, the WD spectra were compared to an appro- 
priate DA model atmosphere. The latter was determined 
iteratively, where we first fitted a high S/N single spec- 
trum with a grid of model atmospheres created by one of us 
(D. Koester, see next section), then used the best fit solu- 
tion to shift the spectra and average them at zero velocity, 
and finally fitted the average again to determine the best 
template. For WD 1743—132 we fitted the single spectrum 
with the grid and determined all parameters simultaneously. 

For all above fits, we multiplied the templates with a 
3 rd degree polynomial to account for the normalization and 
possible variations with wavelength (see §4.5 for details). 
Our best fits gave typical reduced x 2 values of Xrcd.min ~ 1-2, 
2.2 and 1.6 for the WD, the comparison star, and the velocity 
standard, respectively. Best-fit velocities were determined 
by fitting a parabola to the x 2 values to within 60kms~ 1 
of minimum, with uncertainties taken to be the difference 
in velocity over which \ 2 increased by Xied mm (thus effec- 
tively increasing our uncertainties to account for the fact 
that Xred.min did not equal unity). 

For the Keck spectra, we proceeded similarly. Here on 
the red side, we could not use the UVES spectrum due to a 
gap over the Call triplet, and hence we used instead a T e g = 
6000 K, log g = 4.5 dex model by|Zwitter et al.|(|2004|). As we 



trust the absolute wavelength calibration of this observation 
most (being calibrated relative to telluric emission lines), 
we use this estimate of the velocity below to transform all 
velocities to the barycentric reference frame. 

4.2 Radial-velocity orbit and mass ratio 

In Table[l]we list the measured radial velocities for all tar- 
gets, with barycentric corrections applied. For determin- 
ing the orbit, we folded the barycentric velocities using the 
ephemeris in Table[2] and fitted for a circular orbit keeping 
the orbital period and time of accenting node passage fixed 
to the timing values. The fit gave a radial velocity semi- 
amplitude of -ft'obs = 165 ± 7kms~ 1 and a systemic radial 
velocity of 7 = —50 ± 4kms~ 1 with Xrcd — 1-55 for 16 de- 
grees of freedom. 

The radial velocity of the comparison star in the Gemini 
dataset varied as much as 55kms _1 which is considerably 
higher than the uncertainties of individual points. We found 
no evidence for binarity and thus we attribute the large scat- 
ter to systematics, likely induced by slit positioning errors 
and differential atmospheric diffraction. For that reason, we 
chose to use velocities relative to the comparison star, Av. 
This choice relies on the assumption that both the WD and 
the comparison star are affected by the same systematics. 
This should be correct to first order, but given the relatively 
large separation of the two stars on the slit, their different 
distances from the centre of rotation of the instrument, and 
their different colours, small second-order differences may 
remain. Even if any are present however, they should not 
be correlated with orbital phase (since our measurements 
are taken on many different nights), and thus be taken into 
account automatically by our rescaling of the measurement 
errors such that reduced x 2 equals unity. 

After subtracting the velocity of the comparison star, 
we obtain K b B = 166±6kms _1 , A7 = — 101±4kms _1 with 
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Figure 2. Radial velocity measurements of the companion to 
PSR J1738+0333 as a function of the orbital phase. Filled black 
circles depict the points used to fit the orbit and the blue line 
the best-fit solution. Red crosses indicate two outliers that we 
excluded and the black line the best-fit solution with these points 
included. The latter agrees well and is almost indistinguishable. 
The blue triangle shows the Keck point. The green line depicts the 
velocity of the pulsar as inferred from radio timing. All velocities 
are relative to the comparison star, but corrected for its estimated 
61 kms -1 barycentric radial velocity. All error bars represent la 
uncertainties. The orbit is depicted two times for clarity. 



A.rcd 



= 1.07. This orbit is shown in Fig. [2] This fit has two 
outliers, which both are from spectra taken in the night with 
the worst condition (they are also outliers in the relative 
flux between the WD and the comparison star; see Table [T]). 
Excluding these, we find K ^ B = 167 ± 5 kms - and A7 = 
-103 ±3 km s" 1 with xLi = 0.93 for 14 de grees of freedom. 
We will use these latter values as our best estimates, but 
note that all fits gave consistent results, so our inferences do 
not depend on this choice. 

Because the exposure time is a significant fraction of 
the orbit (t oxp ~ 0.12 F\,), the observed semi-amplitude is af- 
fected by velocity smearing. This reduces the measured am- 
plitude by a factor sin(7ri cxp /Pb)/(7rfcxp/-Pb) = 0.976. Thus, 
the true radial- velocity amplitude is T^wd = 171 ±5 kms -1 . 

Likewise, the semi-amplitude of the pulsar's projected 
radial velocity is K PS r = 2ncx/P h = 21.103059(2) kms" 1 , 
where x is the projected semi-major axis of the pulsar orbit. 
Based on the two values calculated above we derive a mass 
ratio of q = K W o/K PS r = 8.1 ± 0.2. 



4.3 Systemic velocity 

The systemic velocity A7 derived above is relative to the 
comparison star. Thus, for an absolute value one needs to 
obtain an estimate of the true velocity of the latter. From the 
Gemini spectra we derived an average value of 64±5kms -1 . 
As discussed above, the individual velocities have a large 



scatter and one may thus worry about systematics. It seems, 
that these are of the order of 15 — 20 kms -1 . First, for 
the velocity standard WD 1743— 132 we find a velocity of 
— 58.6 ± lkms" , which is offs et by 14.2 k m s~ from the 
catalogue value of —72.8 km s -1 (Reid 1996 1. Second, for the 



comparison star, our Keck spectrum yields 61 ± 5 kms - . 
As mentioned above, we believe the wavelength calibration 
is most reliable for the Keck spectrum, so we adapt this 
velocity. For PSR J1738+0333, correcting for the gravita- 
tional red-shift of the white dwarf of 3 kms -1 (using the 
mass and radius derived in §|4]6), we infer a systemic veloc- 
ity of 7 = -42 ± 16 kms -1 . 



4.4 Interstellar reddening 

We calculated the run of reddening along the line of sight us- 



ing the Galactic extinction model of Drimmel et al. ( 2003 1 



We find that the interstellar extinction increases smoothly 
to reach a maximum value of Ay = 0.56 at 1.3 kpc and 
remains constant thereafter. This is similar to the maxi- 
mum value along this line of sight of Ay = 0.65 inferred 
from the maps of Schlegel et al. (19981. Therefore, for 



both PSR J1738+0333 and the comparison value we adopt 
Ay = 0.56 ± 0.09, with the uncertainty taken to be the 
difference between the two models. 

We can now use these results to estimate the distance 
of the comparison star: Adopting My = 4.3 and (B — V)q = 
0.57 for a G OV star ( |Cox|2000[ ) and A B = 1.321,4^ ( |Schlegel 
et al. 1998) we obtain a distance of ~ 4.3 kpc for both bands. 



As a sanity check for the systemic velocity derived above, 
we can calculate the expected velocity of the comparison for 
the photometric parallax: Assuming the Galactic potential 
of |Kenyon et al.| ( |2008[ ), a distance to the Galactic center of 
8.0 kpc and a peculiar velocity of the Sun relative to the lo- 
cal standar d of rest of (U, V, W) = (10.00, 5.25, 7.17) kms" 1 
| [CroE|[2000l ), we find that the local standard of rest at the 
position of the comparison star moves with a speed of 
~ 60 kms -1 . Given the uncertainties of the model and our 
measurements and the possibility of peculiar motion, the 
latter agrees well with our estimated value. 



4.5 Temperature and surface gravity of the WD 

The zero-velocity average spectrum (Fig.[3| shows deep 
Balmer lines up to H12, typical for a WD with a hydrogen 
atmosphere and low surface gravity. 

Quantitative estimates for the atmospheric parameters 
were obtained by modelling the spectrum with a grid of 
DA model atmospheres extending from 7000 K to 25000 K 
and log g = 6.00 to loggr = 8.00 with step-sizes of 100 K 
and 0.1 dex respectively. The models used in this work are a 



recent update of the grid presented in Koester ( 2008 1 which 



incorporates the improved treatment of pressure broadening 



of the absorption lines by Tremblay & Bergeron ( 2009 1 



At each point of the grid that we scanned, we fitted for 
the normalization with a polynomial function of the wave- 
length. This was found necessary in order to account for the 
(up to) ~ 10%, slowly varying continuum deviations, caused 
by in-perfect flux calibration. Assuming our flux calibration 
is perfect (namely, using a normalization factor that does not 
vary with wavelength) resulted in a poor fit with large scale 
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Figure 3. Left: From top to bottom: The zero-velocity flux-calibrated average spectrum of PSR J1738+0333 obtained with Gemini, the 
corresponding best-fit atmospheric model, the (single) spectrum obtained with Keck, and the residuals from the fit (see jj4|. The model 
and the Keck spectrum are shifted down by 10 and 20 units respectively. Right: Details of of the Balmer series in the average spectrum 
(H/3 to H12, from bottom to top), with the best-fit model overdrawn (red lines). Lines are shifted by 8 units with respect to each other. 



structure in the residuals and lines systematically deeper 
than the best-fit model (best-fit values: T e s = 9010 ± 50 K, 



logp = 6.81 ± 0.12 dex with %t 



9). Similarly underesti- 



mated lines were obtained using a fitting routine normally 
used by one of us (D. Koester) that assumes a fixed slope for 
the continuum over the length of each line. The former com- 
parison revealed that there was also a smaller spectral range 
between 4400— 4780 A with features similar with the ones 
seen in the flat fields (see Sj3]| , likely associated with the holo- 
graphic grating (we were alerted to this effect because it was 
much more obvious for the companion of PSR J1909— 3744; 
vK+12). Fortunately, no Balmer lines are present in this re- 
gion, and hence we simply modelled the spectrum excluding 
this range (specifically, we fitted the ranges 3700-4400, and 
4780-4960 A). Like for our radial- velocity fits, we accounted 
for the spectral resolution by convolving the models with a 
truncated Gaussian. 

Using the choices described above we obtain T c g — 
9129 ± 20 K (implying a spectral type DA5.5) and loggr = 
6.55 ± 0.07 dex with Xred mm — 1-5 (f° r ~ 800 points and 
6 parameters). Here, the best-fit values and statistical un- 
certainties were determined by fitting the x 2 surface with a 
paraboloid as in |Bassa et al.| ( |2006| ). We verified these es- 
timates using a Monte-Carlo simulation with 10 6 iterations 
(see Fig.|4|. The results are almost identical, with the sim- 
ulation giving slightly larger uncertainties. However, as we 
will see below, the systematic uncertainties are larger. 

The best-fit model is shown in Fig.[3] Most lines are 
matched almost perfectly, but Hll and H12 are slightly un- 
derestimated. We do not know the reason for this. As the 



continuum matches very well, it cannot be due to errors in 
the flux calibration (which would be multiplicative), while 
most other observational issues (scattered light, etc.) would 
lead to lines that have reduced rather than increased depth. 

Given the above discrepancies, as well as previous ex- 
perience with fitting model atmospheres, it is likely our un- 
certainties are dominated by systematics rather than mea- 
surement errors. We investigated this in three ways. First, 
we tried small changes in the assumed spectral resolution 
(by 5 %) and varied the different polynomial degrees for the 
continuum (2nd to 4th order). The former had only very 
small effect (~20K changes in T and ~ 0.03 dex changes in 
log g ) , while changing the degree of the polynomial caused 
larger differences: 0.1 dex (1.5cr) for the surface gravity and 
up to 150 K (7a) for the temperature. Our central values are 
based on a 3rd degree polynomial, since it gave the best fit 
for the higher lines. 

As a second check, we obtained an independent mea- 
sure of the atmospheric parameters using the Keck spec- 
trum. Again using a third-degree polynomial for the con- 
tinuum, and fitting the same wavelength regions, we find 
T eS = 9281 ± 110 K and log g = 6.57 ± 0.13 dex. Here, 
switching between polynomials for the continuum normal- 
ization had a slightly smaller impact on the estimated val- 
ues (~ 100 K in T and ~ 0.1 in loggr ). While the surface 
gravity agrees almost perfectly with the Gemini value, the 
effective temperature is somewhat higher, suggesting, again, 
that temperature is more sensitive to our modelling assump- 
tions. 

Finally, we fitted the individual spectra with the model 
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Figure 4. Left: Constraints on the temperature and gravity of the white dwarf companion to PSR J1738+0333 inferred from our 
model-atmosphere fit, with contours at A\ 2 = X? c d, m in > 4 X?cd,min> 9x ?cd,min' and 16 X 2 c d,mi[L: Tne horizontal and vertical sub-panels 
show the histograms of the distributions for T e g and log g from our Monte Carlo simulation (9414). Right: Constraints on the mass and 
radius of the WD. The shaded area depicts the distribution of realizations from our Monte-Carlo simulation (fj4~|6). Overdrawn are: the 
central value and ler confidence limits of the observed surface gravity (red lines); the model tracks of |Panei et al.| (|2000|> for constant 
temperature (8000 and 12000 K; dotted lines); the mass-radius relations of Scrcnclli et al. (200^1 (solid) and |Panei et al.| | |2007| l (dashed) 
for our best-fit temperature of 9130 K; errorbars showing the independent constraints from photometry and radio timing. The horizontal 
and vertical panel show the inferred distributions for the WD radius and mass, respectively, as well as the independent photometric 
estimate for the former (in green; see ^4]5). 



atmospheres and obtained a mean temperature of < T >= 
9153 ± 38 K with an rms scatter of 155 K. 

From the above, it is clear the formal uncertainty on 
especially the temperature is too small, and we adopt as 
realistic estimates T e g = 9130 ± 150 K and loggr = 6.55 ± 
0.10 dex. Fortunately, the effect of the larger temperature 
uncertainty on the derived masses is small (see fj4|6), because 
the mass-radius relation is much more sensitive to surface 
gravity than to temperature. For our mass calculation below, 
we thus choose to inflate the original \ 2 map to include 
the systematics mentioned above but preserve information 
about the covariance between parameters. 

Finally we searched the average spectrum for signatures 
of rotational broadening. For that we proceeded in two ways: 
First, we broadened a 9000 K, logy = 6.5 model atmosphere 



radius. In terms of magnitudes, 



using the analytical profile of Gray ( 2005 1 with a limb dark- 



ening coefficient of 0.3 and scanned a grid of rotational veloc- 
ities ^ VrS'mi ^ 1500 km s -1 in steps of 50km s _1 . Second, 
we let all parameters free. In both cases we accounted for 
the spectral resolution of the instrument as above. We find 
the rotational broadening consistent with zero with the la 
upper limit being 440 and 510 km s -1 respectively. 



4.6 White dwarf radius from photometry 

We can use the best-fit atmosphere model, the observed 
fluxes, and the distance to obtain an estimate of the WD 



-51og(d/10pc)-A A = 43.234-5 log(7?/R )-2.5 log Fa+ca 

(1) 

where mx is the apparent magnitude in band A, the nu- 
merical term is — 51og(Ro/10pc), Fa is the emitted flux 
per unit surface area integrated over the relevant fil- 
ter, and ca the zero-point. Convolving the best-fit model 



with the B and V band passes of |Bessell| ( 1990 1 yields 
F B = 6.289 x 10 7 ergcm- 2 s _1 A -1 and F v = 4.353 x 
10 7 erg cm -2 s _1 A -1 . Here the uncertainty due to the fit 
is ~ 5% (mostly due to the ~1.5% uncertainty in tempera- 
ture). Using the zero-points of Bessel (1990), cb = —20.498 
and cv = —21.100, and the reddening infered in ^4]4 we ob- 
tain radii R = 0.042 ± 0.004 R and R = 0.042 ± 0.004 R 
for B and V, respectively (with the uncertainty dominated 
by the uncertainty in the parallax). 



4.7 Masses of the white dwarf and the pulsar 

The mass of the WD can be estimated using a mass-radius 
relation appropriate for low mass helium white dwarfs. We 
use the finite-temperature relation for low-mass WDs from 



Panei et al. ( 2000 1, which gave good agreement for the com- 
panion of PSR J1909-3744 (vK+12). 

For the calculation we proceeded as follows: We sam- 
pled the inflated \ 2 surface derived in i|4]3 in a Monte-Carlo 
simulation using 10 6 points uniformly distributed in the 
left — log g plane. For each point within the expectations, we 
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linearly interpolated the 8000 and 12000 K models of |Panei] 
et al. ( 2000 I for WDs with extended hydrogen envelopes to 



the given temperature and calculated the mass and radius 
at the cross-section of the observed value (which scales as 
g = GM/R 2 ) and the model. Subsequently, we calculated 
the mass of the pulsar, assuming a normal distribution for 
the mass ratio with q = 8.1 ± 0.2 (see ^4]2). Furthermore, 
we calculated the inclination using the mass function /m of 
the binary (sin 3 i = / m (M W d + M PSR ) 2 /M^ D ). 

We show the mass distribution in Fig. [4] Since the mass- 
radius relation is steeper towards higher masses, the com- 
panion's mass distribution is asymmetric, with larger wings 
towards higher masses. The same holds for the distribution 
for the radius, with larger wings towards smaller radii. The 
error on the pulsar mass is dominated by the uncertain- 
ties in the companion's mass estimate. To summarize, the 
values that we will be using for the rest of this paper are: 



fl/f fllSI +0-007, +0.017 y j. j _ -, , 7 +0. 07, +0.14 ,» 

M W D — 0.181_ . 05,_0 .013 M 0> -Mpsr — lAl -0.06,-0.08 M ©> 



7?wd = 0.037 



+0.004, +0.007 
0.003,-0.006 



R0 and i = 32.6°^ 1 °lt, 2 1 1 . Here, 
the errors separated by commas are the corresponding 68% 
and 95% intervals spanned by the Monte-Carlo realizations. 

Finally, we also derived mass estimates using two 
different sets of tracks, that gave reliable results for 



PSRJ1909-3744 (vK+12): The tracks of |Serenelli etal 
(200l| yielded M W d = 0.183+° °°^° °^ M and R W v = 



+0. 005, +0.007 
.00. 

of 



0.037 

tracks 

ent values: Mwd 



Ftp 



almost identical to the above. The 
Panei et al.| (|2007| yielded slightly differ- 

M© and Rwo = 



0.175 



+U.UJ7+, 0.029 
-0.005,-0.006 - 



0.038 



+0.005, +0.010 



R©. However, we note that these models 



-0.003,-0.004 

predict a cooling age much smaller than the characteristic 
age of the pulsar (see next section). 



4.8 Cooling age 

We compared the absolute photometric magnitudes in B and 



V with the theoretical cooling tracks of Serenelli et al. ( 2001 1 



for solar metallicity progenitors to infer the cooling age of 
the WD. We did this by minimizing a % 2 merit function 
based on the sum of differences between observed and model 



fluxes in both bands. The track of Serenelli et al. ( 2001 1 clos- 



est in mass to the companion of PSR J1738+0333 is that 
of a 0.169 M©, for which we find r c ~ 4.2 Gyr. For that 
age and mass, the predicted temperature and surface grav- 
ity are T e ff ~ 8500 K and logy ~ 6.35 dex. For our best-fit 
spectroscopic estimates the same track yields r c ~ 2.6 Gyr. 
Since the observed mass is slightly heavier, its cooling age 
must be somewhat lower. Using the 0.193 M© track, we 
get r c ~ 600Myrs. The large difference is due to the di- 
chotomy around 0.2 M© expected between WDs with thick 
and thin hydrogen atmospheres. Using the tracks of |Panei] 
et al. (20071, for the mass of 0.175 M© inferred using those, 
we again find short ages, r c ~ 500 Myrs from the photometry 
and r c ~ 450 Myrs for the spectroscopic parameters. 

Finally, the suggested relatively large age of the system 
(4 Gyr plus 2-10 Gyr for the progenitor to have evolved ) mo- 
tivated us to compare our observations with models for lower 
metallicity progenitors. Using the 0.183 M©, Z = 0.001 track 

5 Gyr. 



of Serenelli et al. ( 2002 I we obtain r, 



The above analysis demonstrates that with the current 
set of observations it is difficult to constrain the cooling age 
of the WD, since this depends on both the thickness of the 



WD envelope and the metallicity of its progenitor. Future, 
more precise constraints on the parallax and consequently 
on the radius, might help to discriminate between different 
cases. 



4.9 3D velocity and Galactic motion 

In §[2] we computed the two components of the transverse 
velocity based on the parallax and proper motion estimates 
from radio timing measurements of the pulsar. Combined 
with the systemic radial velocity 7 = —42 ± 16kms _1 from 
the optical observations of the white dwarf (Table [I]), we 
have the full 3D velocity and can compute the Galactic 
path back in time (like was done for PSR J1012+5307 by 
|Lazaridis et al.||2009[ ). For our calculations we have used 



the Galactic potential of Kenyon et al 



2008 ) , verifying our 



( 1989 \ and Paczyn- 



results with those of Kuijken & Gilmore 
[ski] ( fl990t l. We infer that the PSR J1738+0333 system has 
an eccentric orbit with a Galacto-centric distance between 6 
and llkpc, and an oscillating Z-motion with an amplitude 
of lkpc and a (averaged) period of 125 Myr. We also cal- 
culated the peculiar velocity of the system with respect to 
the local standard of rest at every transition of the Galactic 
plane (Z = 0) during the last 4 Gyr, and find that it ranges 
between 70 and 160kms _1 . We will discuss this further in 

!Hli- 



5 RAMIFICATIONS 

In Table [2] we list the properties of the system derived in 
previous sections and in Fig. [5] we show our constrains on 
the masses. In what follows we discuss the ramifications of 
our work for stellar and binary astrophysics. 



5.1 Kinematics 

PSRJ1738+0333 has a velocity of 85 ± 17kms _1 with re- 
spect to the local standard of rest that co-rotates with the 
Galaxy (Z = 0) at the distance of the pulsar. The latter com- 
pares well with the mean transverse velocity for the bulk of 
MSPs with measured proper motions (e.g. ~ SSkms" 1 ac- 



cording to Hobbs et al. 2005 1. Our semi-quantitative analysis 



in i|4]8 shows that the system's velocity varies as much as 
150kms _1 over the course of its Galactic orbit. Based on 



the simplified potential of Kenyon et al. ( 2008 1 used herein, 
PSR J 1738+0333 has a peculiar velocity between 70kms~ 1 
and 160 kms -1 when it crosses the Galactic plane (Z = 0). 
Thus, assuming that the system had a small peculiar mo- 
tion before the SN explosion, the systemic velocity after the 
formation of the NS must have been in that range. This is 
consistent with a SN explosion with a small, or even negli- 
gible kick ( |Tauris fc Bailes|1996| |Nice fc Taylor|1995[ ). 



5.2 Evolutionary history 

Millisecond pulsars with low-mass helium WD companions 
are expected to form through mainly two different channels 
depending on the initial separation of the progenitor binary 
(e.g. [Tauris 201l] and references therein). The initial sepa- 
ration of the progenitor binary determines the evolutionary 
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Table 2. Properties of the PSR J1738+0333 system. Numbers 
in parentheses (where given) are the formal errors of the best-fit 
model. For details on the timing analysis, including uncertainties 
not relevant here, see Paper II. 

Timing parameter Value 



Reference Epoch (MJD) 
Time of ascending node (MJD) 
Right ascension, a f ,72000) 
Declination, 6 (J2000) 
fj, a (masyr -1 ) 
W (masyr -1 ) 
Parallax, 7r (mas) 
P (ms) 
P (as- 1 ) 

Dispersion measure, DM (cm -3 pc) 
Orbital period, (days) 
Projected semi-major axis, x (lt-s) 
Eccentricity, e 
Mass function, /(Mq) 



54600 

53400.297958820(6) 
17 h 38 m 53?.965 
03deg33'10'.'866 
+7.058(5) 
+5.176(10) 
0.67(4) 
5.85 

2.412xl0- 20 
33.77 
0.35479 
0.3434 
3.5(1.1) x 10~ 7 
0.0003455012(12) 



Optical parameter 



Value 



Temperature (K) 9130(150) 

Surface gravity (logy, spectroscopy) 6.55(10) 

Surface gravity (logy, + Q + 7r+photometry) . . 6.45(7) 

Photometry, V-band 21.30(5) 

Photometry, B-band 21.71(4) 

Semi-amplitude of radial velocity, i^wD (kms -1 ) 171(5) 

Systemic radial velocity, 7 (kms -1 ) —42(16) 

Transverse velocity, Vx (kms -1 ) 59(6) 

3D velocity amplitude (kms -1 ) 72(17) 

Mass ratio, q 8.1(2) 

WD mass, M WD (M Q , spectroscopy) °' 181 to 005 

WD mass, M W D (M©, q + P h ) 0.182 ± 0.016 

WD radius (Spectroscopy) (R ) 0.037±[J;™3 

WD radius (Photometry) (R G ) 0.042(4) 

Cooling age, t c (Gyr) 0.5 - 5 

Pulsar mass, A/ PS r (M g ) l^t™ 7 ; 

inclination, i (degrees) 32.6(1.0) 



status of the donor star at the onset of the Roche lobe over- 
flow (RLO): 

• Case A RLO: For systems with initial periods short 
enough to initiate mass transfer on the main-sequence, it 
is expected that magnetic-braking (aided to some extent by 
gravitational radiation) drives the system to shorter periods, 
resulting in a compact binary in an orbit which is close to 
being perfectly circular (the eccentricity, e < 10 -5 ). These 
systems were first studied in detail by Pylyscr fc Savonije| 
(19891. 



• CaseB RLO: For progenitors with larger initial separa- 
tions the mass transfer is expected to start at a later phase, 
since the star fills its Roche lobe only during shell hydro- 
gen burning, while moving-up the red giant branch. In this 
case the orbit will diverge resulting in a wider binary. In- 
terestingly, for systems following this path, there are two 
theoretical predictions that can be verified observationally: 
The first is a correlation between the orbital period and the 
mass of the WD companion which results from the unique 
relation between the radius of the giant donor and the mass 
of its core which eventually forms the WD ( Savonije|1987 1. 
The second is a correlation between the orbital eccentricity 



and the orbital period ( Phinney 1992 1 arising because the 



turbulent density fluctuations in the convective envelope — 
of which the size increases in more evolved stars (wider or- 
bits) — do not allow for a perfect tidal circularization. 



The critical period that separates diverging from con- 
verging systems (often called bifurcation period) is expected 
to be ~ 1 day, however its precise value depends on the treat- 



ment of tidal interactions and magnetic braking (e.g. Pylyser 
fc Savonije|[l989| |van der Sluys et al.||2005} |Ma fc Li||2009 l 



and is still a subject of debate. The residual eccentricity in 
binaries with radiative donors (i.e. those binaries that evolve 
to tight converging systems) should be closer to zero com- 
pared to binaries in wider orbits but it is difficult to estimate 



by how much, as pointed out by Phinney & Kulkarni ( 1994 1 



With a current orbital period of 8.5 h, PSR J1738+0333 
is most likely the fossil of the former case (Case A RLO). 
However it is interesting to note that our mass estimate 
and the non-zero eccentricity derived in Paper II pass both 
tests for the latter case mentioned above (Case B RLO) 
that predict m c = 0.18 ± 0.01 M ( |Tauris fc Savonije||1999 \ 
and e ~ 4 x 10 -7 , respectively (deduced by extrapolating 



the Phinney 1992 relation to the observed period). This 
apparent agreement seems to be confirmed not only for 
PSR J1738+0333 but also for the other short-period LMWD 



binaries with measured masses (PSR J1012+5307, Lazaridis 
|et al .||2009| PSR J0751+1807, |Nice et al.||2008|), as well as 



low-mass WD companions to non-degenerate stars (e.g. 



Kerkwijk et al. 2010 Breton et al. 20121. Since compan- 



ion masses in converging systems are not expected to fol- 
low these relations, we cannot exclude a coincidence, but 
the matches seem to suggest that there is a grey zone with 
properties from both cases - something which should help 
improve our as yet rather simplified models of these systems. 



5.3 Pulsar mass and efficiency of the mass 
transfer 

Regardless of the evolutionary path followed, the mass trans- 



fer was sub-Eddington (e.g. Tauris & Savonije 1999 1 and 
thus one would expect that a substantial fraction of the 
mass leaving the donor was accreted by the neutron star. 
For PSRJ1738+0333, this is demonstrably false: The min- 
imum mass of the donor star can be constrained from our 
WD mass estimate to be > 1 M a because the available time 
for evolution is limited by the Hubble time (minus the cool- 
ing age of the WD). The amount of mass lost by the donor 
is Mdonor — A^wd, while the amount accreted by the pulsar 
is Mpsr — Afpsiij with the last term being the birth mass 
prior to accretion. For any realistic birth mass of the neu- 
tron star at the low end of its "canonical" birth mass range 
(> 1.20 M Q ), we find that more than 60% of the in-falling 
matter must have escaped the system (after correcting for 
the conversion from baryonic mass to gravitational mass). 
This translates to an accretion efficiency of only e < 0.40. 



This result confirms the findings of Tauris & Savonije ( 1999 1 



who concluded that a substantial fraction of the transferred 
matter in LMXBs is lost from the system, even at sub- 
Eddington mass-transfer rates. 

More constraining (but less stringent) estimates are 
also obtained for the 6.3 h orbital period binary, 



PSRJ0751+1807 (Nice et al.||2008) for which we find s ~ 
0.1 - 0.3. 

Possible mechanisms for mass ejection discussed in the 
literature include propeller effects, accretion disc instabili- 
ties and direct irradiation of the donor's atmosphere from 
the pulsar (e.g., Illarionov & Sunyaev 1975 van Paradijs 
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Figure 5. Constrains on the masses derived from our observa- 
tions. The shaded region is excluded because it would require 
sinj > 1. Our constraint on the mass ratio is indicated by diago- 
nal lines, and that on the WD mass by horizontal lines, with the 
red ones giving the central value and la uncertainties using the 
mass-radius relation of Panci ct al. (2 000[ |, the black dashed ones 
the la range inferred using the models of|Panei et al,||2007|l and 
the black dotted the la limit based on |Serenelli et al.| ( |2001[ ). 



1996 Dubus et al. 19991. Alternatively, the neutron star 



in PSR J1738+0333 might have formed via the accretion- 



induced collapse of a massive ONeMg WD (e.g., Nomoto & 



Kondo 19911. If the neutron star was formed towards the 



end of the mass transfer it would not have accreted much 
since its birth. A possible problem with the above mecha- 
nism however, is that it is specific to pulsars, while similarly 
inefficient accretion has been found also for low-mass WDs 
with non-degenerate companions (e.g., |Breton et aL||2012[ ), 
suggesting the problem in our understanding is more gen- 
eral. Finally, we note that even major inefficiencies in the 
mass accretion process do not pose a problem for the recy- 
cling scenario: the accreted mass needed to spin-up a pulsar 
to a ~ 5 ms period is only of the order of 0.05 Mq. 



6 CONCLUSIONS 

The main result of this work is the determination of the 
component masses of the PSR J1738+0333 system and adds 
to the three previously known MSP-LMWD binaries with 



spectroscopic information (PSR J1012+5307, van Kerkwijk 
etaLl(p96l) and|Callanan et al.| ( fl998| ); PSR J1911— 5958A, 
Bassa et al.| ( |2006[ ); PSR J1909— 3744, vK+12). 

Our mass estimates are derived independently of any 
strong field effects and thus transform the PSR J 1738+0333 
system into a gravitational laboratory, which - due to its 
short orbital period, gravitationally asymmetric nature, and 
timing stability - provides the opportunity to test the ra- 
diative properties of a wide range of alternatives to GR (see 
Paper II for details). 



Based on our measurements of the component masses, 
GR predicts an orbital decay of A = -2.77^°;^ x 10" 14 . 
While the actual P\, inferred observationally is still less pre- 
cise than this prediction, it will eventually provide a precise 
test for the input physics of atmospheric and evolutionary 
models. Assuming the validity of GR, one can confront the 
spectroscopic WD mass estimate implied by the mass ra- 
tio and intrinsic orbital decay of the system and thus test 
the assumptions for stellar astrophysics and WD composi- 
tion that were used to model the evolution of the WD. Ad- 
ditionally, this mass estimate, combined with parallax and 
absolute photometry constrains independently the surface 
gravity of the WD. The current estimates on these parame- 
ters imply a surface gravity of \ogg = 6.45 ± 0.07 dex. While 
this is formally more accurate than our spectroscopic con- 
straint, it might still be dominated by systematics on the 
distance, arising from correlations between the parallax and 
DM variations (see Paper II for details). 

Finally, the interpretation of the mass estimates within 
the context of our current understanding for binary evolu- 
tion implies that a significant fraction of the accreted mate- 
rial during the LMXB phase is ejected by the system. The 
discovery and study of more similar systems in the future 
will allow further tests of this result. 
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